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In this talk we present a model to demonstrate how time-periodic potential can be used to ma- 
nipulate quantum metastability of a system. We study metastability of a particle trapped in a well 
with a time-periodically oscillating barrier in the Floquet formalism. It is shown that the oscillat- 
ing barrier causes the system to decay faster in general. However, avoided crossings of metastable 
states can occur with the less stable states crossing over to the more stable ones. If in the static 
well there exists a bound state, then it is possible to stabilize a metastable state by adiabatically 
increasing the oscillating frequency of the barrier so that the unstable state eventually cross-over to 
the stable bound state. It is also found that increasing the amplitude of the oscillating field may 
change a direct crossing of states into an avoided one. Hence, one can manipulate the stability of 
different states in a quantum potential by a combination of adiabatic changes of the frequency and 
the amplitude of the oscillating barrier. 

I. INTRODUCTION 

Ever since the advent of quantum mechanics, quantum tunneling has been an important and fascinating subject. 
This phenomenon arises frequently in physics. In fact, one of the first successful applications of quantum mechanics has 
been the explanation of the a-decay of atoms as a quantum tunneling process [l| . Recent examples include tunneling 
phenomena in semiconductors and superconductors @, in Josephson junction systems resonant tunneling in 
hetero junction nanostructures [1], tunneling ionization of atoms photon-assisted tunneling in superconducting 
junctions and semiconductor superlattices [fj, etc. 

In cosmology, quantum metastable states play an essential role in some versions of the inflationary models of the 
early universe Q . In these models inflation of the early universe is governed by a Higgs field trapped in a metastable 
state. Inflation ends when the metastable state decays to the true ground state of the universe. During inflation 
the universe expands exponentially. It is thus obvious that the metastable state of the Higgs field is trapped in a 
rapidly varying potential. The problem is therefore a truly time-dependent one. Unfortunately, owing to the inherent 
difficulties of the problem, more often than not one has to consider the decay of the Higgs field in a quasi-stationary 
approximation, in which the decay is studied assuming a static potential (8j. Before one can deal with the decay 
of Higgs field in the non-stationary potential, it is desirable to gain some insights first by studying metastability in 
time-dependent potential in simple quantum-mechanical models. 

An early attempt at studying the effects of time- varying forces on quantum metastability appears in Fisher's work 
Q, which was motivated by an experiment on quantum tunneling of the phase in a current-biased Josephson junction 
with a weak microwave perturbation [10]. In this work Fisher considered the general problem of quantum tunneling 
in a metastable well with a weak oscillatory force. There he reformulated the standard WKBJ approach to quantum 
decay in order to include a weak time-dependent perturbation. For a class of metastable potentials which interpolates 
between the cubic potential and a truncated harmonic-oscillator potential, he showed that the decay rate is generally 
enhanced by the weak oscillatory force. The potential considered by Fisher has a number of oscillator-like levels near 
its minimum. The opposite situation where only two levels are present was considered by Sokolovski [Tl| , who studied 
the effect of a small AC field mixing two levels in the well on the tunneling rate in a semiclassical framework. 

The results in 0, [ll[ are quite general for a class of weak oscillatory forces. However, it is desirable to consider 
other possibilities, e.g., exact solutions and/or nonperturbative results. 

In [l2l . [l3j we have considered simple driven quantum metastable models in which a particle is trapped in a well 
with a periodically driven rectangular barrier. In order to do away with any restriction of amplitude or frequency of 
the driving force, and of the number of states in the potential, we treat the problems in the framework of the Floquet 
formalism [3, 12 EH . An exact expression determining the Floquet quasi-energies of stable/metastable states in the 
well is derived. From the solution of this equation we find that while the oscillating barrier makes the system decay 
faster in general, there is the possibility that avoided crossings of metastable states can occur with the less stable 
states crossing over to the more stable ones. 

That an oscillating potential can affect the tunneling property of a system has also been noticed before, eg. in 
quantum transport process [l||[l7j])- In [n| it was demonstrated that an propagating particle at appropriate incident 
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energy can be trapped into a bound state by an oscillating square-well. In [17fl it was found that a particle can be 
localized in one side of a time-dependent double well if the amplitude and the frequency of the oscillating field were 
chosen properly. Our examples shows how a time-periodic field can modify the metastability of a decaying state. 



II. THE MODEL 



To be concrete, in this talk we present the model considered in [12|. This model consists of a particle of mass m 
trapped in a square well with a harmonically oscillating barrier, 

x < , 

V(x,t) = \ ?' . . °f *<?' (1) 




\ cos(wi) , a < x <b 
x > b . 

Here Vq, Vi, Vq and u> are positive parameters, with V ' < Vq and V\ < Vq — Vq. According to the Floquet theorem, the 
wave function of a time-periodic system has the form v E' £ (x, t) = e~ let / h <& e {x, t), where $ £ (x, t) is a periodic function 
with the period T = 2it/uj, i.e. Q e (x, t + T) = $ e (x, t), and e is the Floquet quasi-energy, which we will call Floquet 
energy for brevity. It should be noted that the Floquet energy is determined only modulo Uu>. For if {e, $ e } is a 
solution of the Schrodinger equation, then {e 1 = e + nhu , $ e / = $ e exp(inujt)} is also a solution for any integer n. 
But they are physically equivalent as the total wave function \l/ e is the same [15| . All physically inequivalent states 
can be characterized by their reduced Floquet energies in a zone with a width Tim. We therefore consider solutions of 
e only in the first Floquet zone, i.e. e 6 [0, Two). 

The wave function of the particle is described as follows [l2j : 

= e- ist / n $ e (x,t) (2) 

f E"=-oo A. sin(fc„a ; )e- m " t , < x < a , 

= e- iEt / fi <^ E^-oo E~-oc (W qix + he-* x ) J n -i (t) e ' mut > <* < » < b , 

where 



fc„ = y/2m(e + nhu))/h , 

qi = y/2m(V - e - lhu)/h , (3) 



k' n = ^2m(£ + nhuj -Vq)/% , 



and J„'s are the Bessel functions. In the region x > 6, we have adopted Gamow's outgoing boundary condition, 
namely, there is no particle approaching the barrier from the right Matching the wave function and its first 
derivative at the boundaries x = a and x = b, we obtain the relations among the coefficients A n , a n , b n and t n : 

A n sm(k n a) - Y,( a i eq ' a + b i e ~ q ' a ) J »-'(«) ' 
i 

k n A n cos(k n a) = J^H ( a ' e?!a ~ bte- q ' a ) J„_ ; (a) , (4) 
i 

t n e lk > - Y,( a ^ q ' b + b ^ qib ) J «-'( a ) > 
I 

ik' n t n e lk '~ b =Y,qi (aie<>' b - b ie ^' b ) J„_ ; (a) , 
l 

where a = Vi/huj. The Floquet energy is determined from these relations by demanding non-trivial solutions of the 
coefficients. In practice, however, we must truncate the above equations to a finite number of terms, or sidebands as 
are usually called in the literature, eg. n = 0, ±1, . . . , ±N . The number N is determined by the frequency and the 
strength of the oscillation as N > Vi/Tiuj 

By demanding non-trivial solutions of the coefficients a n , b n , A n and t n in Eq.([4]), we obtain, after some tedious 
algebra, an equation which determines the Floquet energy e: 

T-i Qo i , , jp F$q Q +iF 6 k' ( qo \ -2„ (b-a) fr\ 

F i —tank a a + F 2 = — — — [ F 3 — tanfc a - F x e qo{ >. (5) 
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Here Fi(ko, k' Q ,co, Vi) are functions of fco, <Zo, k' (see [l2| for their complete forms), and hence are connected with the 
Floquet energy e (c.f. Eq. ([3])). If the solutions e of Eq.© are complex (real) numbers, the corresponding Floquet 
states are metastable (stable) states. The non-decay probability P(i), which is the probability of the particle still 
being trapped by the potential barrier at time t > 0, is given by 



P(t) 



foMx,t)\ 2 dx 

i r. I i 1 2 tl r 

(6) 



/ 6 |*(x,0)| 2 ^ 
r 2i m ( E )t/n Jo \M^,t)\ 2 dx 



1Im{e)t 



/h h(t) , 



with -P(O) = 1. The imaginary part of the Floquet energy, which enters P(t) via the factor exp(2Im(e)t /%) , gives a 
measure of the stability of the system. Unlike the static case, however, here P(t) is not a monotonic function of time, 
owing to the time-dependent function h(t) after the exponential factor in Eq.®. But since h(t) is only a periodic 
function oscillating between two values which are of order one, the essential behavior of P(t) at large times is still 
mainly governed by the exponential factor. Hence, as a useful measure of the non-decay rate of the particle in the 
well, we propose a coarse-grained non-decay probability P(t) defined as 

Pit) =e 2Im ^ h {h(t)) , (7) 

where (h(t)) is the time-average of h(t) over one period of oscillation (some graphs of P(t) are presented in [l2l|). 

The coefficients Fi(ko, k' , lo, V\) all approach to one in the limit a = Vi/Tiio — > 0. Hence in the limit V\ — ► or 
lo — » oo, Eq.© reduces to the corresponding equation for the case of static potential with potential Vo in the region 
a < x < b, and the Floquet energy in this limit is just the (real or complex) eigenenergy of the static case. This is 
understandable, since in the limit V\ —* the potential becomes static, and at high frequencies the particle in the 
well will only see a time- averaged barrier of effective height Vq. 



III. NUMERICAL RESULTS 



We now study numerical solutions of Eq.© with a specific potential. We take a = 1, b = 2, U = 15 and V ' = Vq/2 
in the atomic units (a.u.) (e = m e = h = 1). In the static case this potential supports one bound state, with energy 
E Q /V Q = 0.232123, and one metastable state, with complex energy E x /Vq = 0. 864945-0. 002552611 For the oscillating 
potential, we solve Eq.© in 2-sideband approximation, i.e. we take N = 2. This is accurate enough for oscillating 
frequency to > Vi/2. 

Fig. 1 and 2 present the graphs of the real and imaginary parts of the Floquet energy (e/Vb) as a function of 
w/Uo > 0.2 with Vi = O.lVb and 0.2Vo, respectively. We find that the solutions of Eq.© have the form e = Eo + nuj 
(n = 0, ±1, ±2, . . .), with K(eo) (the horizontal branch) lie close to the energies Eo and 3?(i?i) in the static potential. 
That is, these branches of 5ft(e) emanate from either Eq or K(£'i) at u> = 0. Branches emerging from the same point 
have the same imaginary part. Numerical results show that, with the barrier oscillating, the stable state (Eq) in the 
static case becomes unstable, and the unstable state (Ei) will decay even faster. For simplicity, in Fig. (la) and (2a) 
we show only six branches (n = 0, ±1, ±2 and —3) emerging from K(£'i), and only the central branch (n = 0) and a 
subband (n = —1) from Eo. As mentioned before, we only take solutions in the first Floquet zone, 5Ff(e) (modulo uS), 
which are points under the line 5R(e) = cj. 

From Fig. 1 and 2 we also see that a direct crossing occurs at frequency lo ~ 5R(i?i — Eo)/2 (point c). However, 
independent of the values of V\, as u) approaches the frequency u ~ 5R(£'i — E ) = 0.632822Vb, an avoided crossing 
(e, e') between the real parts of the Floquet energies occurs. Fig. 2 indicates that larger values of V\ only enhance 
the instability of the system and the repulsion between the two levels at avoided crossing. Thus as the frequency lo is 
increased, the state emanating from E\ has Floquet energy with real part given by values along the path abb'cdd'e' f 
(the dark dotted curve), while the real part of Floquet energy of the state emerging from Eo lies along the path 
cegg'h (the solid curve). The imaginary parts of these two paths are depicted in Fig. (lb) and (2b). One sees that an 
exchange of the imaginary parts takes place at the avoided crossing ee'. Beyond the avoided crossing, the upper state 
becomes more stable than the lower state. This gives the possibility of stabilizing an unstable state by an oscillating 
field. We recall that as lo — > oo, Eq.© reduces to the one for the static potential. In the example considered here, 
the lower state supported by the static well is a stable bound state, and hence the unstable upper Floquet state can 
be made stable in the high frequency limit. Even more simply, the same aim can be achieved by adiabatically tuning 
down the amplitude V\ just after the avoided crossing, as in this limit the potential becomes the static one. 
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Finally, it is interesting to note here the role of amplitude V\ of the oscillating barrier in the model (a detailed 
study is presented in [l3| . As we have seen, the presence of a non-vanishing V\ always makes the system less stable. 
However, if V% is reduced, an avoided crossing may turn into a direct one. In the present case, the avoided crossing 
ee' changes into a direct crossing for Vi/Vq < 0.03. Conversely, increasing V\ could change a direct crossing into an 
avoided crossing. At an avoided crossing, the imaginary parts of the Floquet energies cross, while the real parts do 
not. At a direct crossing, it is the real parts, not the imaginary parts, that cross. But the more (less) stable state 
has the tendency to become less (more) stable. This is evident from the Floquet energy at the direct crossing point 
c in Fig. 1 and 2. These observations are consistent with the semiclassical results obtained in [llT | by perturbative 
methods. Hence, by a combination of adiabatic changes of the frequency and the amplitude of the oscillating barrier, 
one can manipulate the stability of different states in a quantum potential: tune up V\ until a direct crossing becomes 
an avoided one, increase u so that the avoided crossing is passed, then reduce V\ to make the potential static. In the 
process, two states in the well are interchanged. 



To summarize, our results show that an oscillating potential barrier generally makes a metastable system decay 
faster. However, the existence of avoided crossings of metastable states can switch a less stable state to a more stable 
one. If in the static well there exists a bound state, then it is possible to stabilize a metastable state by adiabatically 
changing the oscillating frequency and amplitude of the barrier so that the unstable state eventually cross-over to the 
stable bound state. Thus an time-dependent potential can be used to control the stability of a particle trapped in a 
well. 

The model we present here has its barrier oscillating. One may also consider the same model with an oscillating 
bottom instead. However, this latter model turns out to have identical Floquet energy spectrum as the one presented 
here for the same frequency and amplitude of the oscillating field. This can be explained by a discrete transform 
which connects the equations of boun dary conditions of the two models, and by a gauge transformation which maps 
the wave functions of the two models 1131. 




IV. SUMMARY 
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Figure Captions 

Figure 1. The Floquet energies (e/Vb) of the two metastable states versus the barrier oscillating frequency (uj/Vq) 
for Vb = lba.u.,V^ = Vo/2 and Vi/Vo = 0.1 in the atomic units (a.u.) (e = m e = K = 1) . In (a) the real 
parts of the Floquet energies are shown in the first Floquet zone under the line 5R(e) = lo (the straight line). 
The light dotted lines show how the different branches emanate from the two states in the static case (with 
Eq/Vq = 0.232123, and %t(Ei)/Va — 0.864945). In (b) the corresponding imaginary parts of the Floquet energies 
of the two states are plotted. The dotted curve corresponds to the state with real parts given along the path 
abb'cdd'e' f, and the solid curve corresponds to the state with real parts given along cegg'h. 



Figure 2. Same plot as Fig. 1 for V = lha.u., V ' = V /2 and Vi/V = 0.2. 
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